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ABSTRACT 

Tidal debris around galaxies can yield important clues on their evolution. We have identified tidal debris in 1 1 
early type galaxies (T < 0) from a sample of 65 early types drawn from the Spitzer Survey of Stellar Structure in 
Galaxies (S 4 G). The tidal debris includes features such as shells, ripples and tidal tails. A variety of techniques, 
including two-dimensional decomposition of galactic structures, was used to quantify the residual tidal features. 
The tidal debris contributes ^3 - 10% to the total 3.6 fim luminosity of the host galaxy. Structural parameters 
of the galaxies were estimated using two-dimensional profile fitting. We investigate the locations of galaxies 
with tidal debris in the Fundamental Plane and Kormendy relation. We find that galaxies with tidal debris 
lie within the scatter of early type galaxies without tidal features. Assuming that the tidal debris is indicative 
of recent gravitational interaction or merger, this suggests that these galaxies have either undergone minor 
merging events so that the overall structural properties of the galaxies are not significantly altered, or they have 
undergone a major merging events but already have experienced sufficient relaxation and phase-mixing so that 
their structural properties become similar to those of the non-interacting early type galaxies. 

Subject headings: galaxies: structure — galaxies: elliptical and lenticular, cD — galaxies: interactions — 
galaxies: evolution — galaxies: fundamental parameters — galaxies: individual (NGC 474, 
NGC 680, NGC 1222, NGC 1482, NGC 2634, NGC 3032, NGC 4106, NGC 4694, NGC 
5018, NGC 5574, NGC 5576) 



1. INTRODUCTION 
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Early type galaxies (ETGs) used to be regarded as relatively 
simple, quiescent stellar systems undergoing passive evolu- 
tion with little or no star formation. Various studies, however, 
indicate that the histories of ETGs may be more complex. Fol- 
lowing the suggestion that most field elliptical galaxies may 
be the result of merging of disk galaxies (Toomre 1977), nu- 
merous observational studies have shown that many ETGs 
have faint features or disturbances around them that are sug- 
gestive of recent interactions. In particular, many ETGs show 
tidal features such as "shells" (Malin & Carter 1980), "rip- 
ples" (Schweizer 1980) or distinctive tidal tails (Arp 1966; 
Malin & Carter 1983; Fort et al. 1986; Schweizer & Seitzer 
1992; Michard & Prugniel 2004; van Dokkum 2005; Canal- 
izo et al. 2007; Bennert et al. 2008; Tal et al. 2009; Janowiecki 
et al. 2010). Although some shells are slightly bluer or red- 
der than the underlying galaxies, most shell colors are found 
to be similar to those of their underlying galaxies (Fort et al. 
1986; Forbes et al. 1995b; Turnbull et al. 1999; Sikkema et al. 
2007). On the other hand, tidally disturbed ETGs, are on av- 
erage bluer than non-interacting ETGs and the more tidally 
disturbed, the bluer (Schweizer & Seitzer 1992; Michard & 
Prugniel 2004; Tal et al. 2009). 

There have been various studies on the formation of shells. 
Numerical simulations of galaxy mergers of unequal mass, 
mass ratio of 1/10 - 1/100 (i.e minor mergers) produce shells 
(Quinn 1984, Dupraz & Combes 1986, Hernquist & Quinn 
1987a, Hernquist & Quinn 1987b). In these simulations, 
shells are density waves from infalling stars of a low mass 
and low velocity dispersion companion galaxy during a merg- 
ing event. Major merger models (mergers of comparable mass 
galaxies) can also form shells (Hernquist & Spergel 1992, Hi- 
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bbard & Mihos 1995). In these cases, the outer disk material 
of the pre-merger spiral falls into the merged remnant late, 
long after the inner region of the remnant has relaxed in the 
merging event (Hernquist & Spergel 1992). A weak inter- 
action model (Thomson & Wright 1990, Thomson 1991) is 
also able to produce long lasting (^10 Gyr), azimuthally dis- 
tributed shell structures by the interference of density waves. 
In this scenario, the shell material does not originate from 
a companion galaxy, but comes from the thick disk popu- 
lation of dynamically cold stars of the galaxy. This is an 
unexpected result because galaxies that have a disk compo- 
nent should have been classified as SOs, rather than elliptical. 
Thus, it was difficult to explain the formation of shells in el- 
liptical galaxies with a weak interaction model. But using 
two-dimensional stellar kinematics, SAURON (Bacon et al. 
2001) and ATLAS 3D (Cappellari et al. 2011a) studies have 
found that ETGs can be divided into two kinematically dis- 
tinct families, slow and fast rotators (Emsellem et al. 2007; 
Emsellem et al. 2011). Fast rotators show large-scale rota- 
tion, while slow rotators show regular velocity field without 
large-scale rotation. These classifications are nearly insensi- 
tive to viewing angles. Emsellem et al. (201 1) show that 66% 
(45/68) of elliptical galaxies are fast rotators, which implies 
that they might have a disk component even though they were 
classified as ellipticals. Thus, at least in this fraction of ellip- 
tical galaxies, shells might be formed via weak interaction as 
well as via major or minor mergers. 

Van Dokkum (2005) reports that as many as 71% (61/86) of 
bulge-dominated ETGs in the NOAO Deep Wide-Field Sur- 
vey at (z) ^0.1 show morphological evidence of tidal inter- 
actions. A similar result was obtained by Tal et al. (2009), 
who found that 73% (40/55) of nearby (15 - 50 Mpc) ellipti- 
cal galaxies show such features. With deep V-band imaging 
data, Janowiecki et al. (2010) investigated the five brightest el- 
liptical galaxies in the Virgo cluster and found that except for 
M84, all showed shells and diffuse tidal streams in their outer 
halos, in the range = 26 - 29 mag arcsec" 2 . At slightly 
higher redshifts, (z) - 0.15 - 0.20, Bennert et al. (2008) stud- 
ied the merger remnants of 5 early type QSO host galaxies and 
found tidal tails and shells in 4 of them. These studies suggest 
that many elliptical galaxies have undergone recent mergers 
or interactions, consistent with Toomre's (1977) suggestion. 
What remains less well-known is the amount of stellar mass 
in the tidal debris around ETGs. Prieur (1988) investigated a 
shell galaxy, NGC 3923, using B-band images and found that 
the B-band luminosity of the shell contributes about 5% of the 
galaxy luminosity. Bennert et al. (2008) also investigated the 
luminosity of shells in four early type quasi-stellar object host 
galaxies with HST V-band F606W images and estimated that 
shells comprises about 5 - 10% of the total V-band luminosity 
of each galaxy. The current study aims to measure the stellar 
mass fraction of the tidal debris using mid-infrared data at 3.6 
\im from the Spitzer Survey of Stellar Structure in Galaxies 
(S 4 G, Shethetal. 2010). 

A second aim of this paper is to understand whether the 
recent merger or interaction, as revealed by the presence of 
tidal debris, has affected the overall structural properties of 
the galaxy itself. Kormendy (1977) showed that elliptical 
galaxies and the bulges of SO galaxies both have correlations 
between effective radius and surface brightness (Kormendy 
relation, hereafter KR). In addition, ETGs are known to fol- 
low empirical scaling relations between the effective radius, 
the surface brightness and the central velocity dispersion, re- 



ferred to as the Fundamental Plane (hereafter FP, Dressier 
et al. 1987; Djorgovski & Davis 1987) or the Fundamental 
Manifold (Zaritsky et al. 2006; Zaritsky et al. 2008) if one 
extends FP formalism to lower or higher masses from dwarf 
spheroidal galaxies up to intracluster stellar populations of 
galaxy cluster. But not all ETGs lie on the FP; post starburst 
galaxies (E+A's) are one examples of galaxies that do not 
(Yang et al. 2008). Similar to ETGs, E+A galaxies have large 
bulge to total luminosity ratios, large Sersic indices, high con- 
centration indices and a lack of on-going star formation, but 
have significantly larger asymmetry indices. However, these 
galaxies do not lie on the same FP with ETGs and follow their 
own relation, a tilted FP (Yang et al. 2008). Yang et al. (2008) 
showed that 55% (11/21) of their E+A galaxies have dra- 
matic tidal features, indicative of a recent merger/interaction. 
There might be two possible explanations for the offsets from 
the standard FP: changes in the structure, and/or changes in 
the underlying stellar populations. Through gas-poor ("dry") 
merging events, galaxies can increase their half light radius 
(Nipoti et al. 2003; Boylan-Kolchin et al. 2006; Naab et al. 
2009; van der Wei et al. 2009; Hopkins et al. 2009) and their 
structural properties change and they may lie off from the FP 
of normal ETGs. On the other hand, changes in the stellar 
populations due to star formation via gas-rich ("wet") merg- 
ing events may cause galaxies to have smaller mass to light 
ratios (M/L) and thus deviate from the standard FP. For ex- 
ample, E+A galaxies are found to have 3.8 times smaller M/L 
(Yang et al. 2008) and stand apart from the standard FP. As 
tidal debris and shells around ETGs are indicative of recent 
merger or interaction, there might be changes in structure of 
those ETGs. Therefore we also address whether the ETGs 
with shells and tidal debris lie off the FP and/or KR. 

The paper is organized as follows: In Section 2, we de- 
scribe the data and data analysis techniques. Then, we de- 
scribe galaxies showing shells, tidal debris and the fraction of 
tidally disturbed galaxies in ETGs in Section 3. The luminosi- 
ties of the tidal debris are presented and discussed in Section 
4, and structural properties, KR and FP of tidally disturbed 
ETGs are explored in Section 5. Finally, we summarize our 
results and state our conclusions in Section 6. Throughout 
this paper, we adopt Ho = 73 km sT l Mpc~ l , Cl m = 0.27, and 
£1a = 0.73 when estimating distances. 

2. OBSERVATIONS AND DATA ANALYSIS 
2.1. Data 

S 4 G is one of the ten large Exploratory Science programs 
underway in the Spitzer post-cryogenic mission. The overall 
goal of S 4 G is to image over 2,300 nearby galaxies at 3.6 and 
4.5 \im through new "Warm" observations obtained with the 
Infrared Array Camera (IRAC, Fazio et al. 2004) instrument 
on the Spitzer Space Telescope (Werner et al. 2004). These 
data will add to the existing data of galaxies from the Spitzer 
Heritage Archive. S 4 G sample galaxies have radial velocity 
V r adio < 3,000 km s" 1 corresponding to d < 40 Mpc, a to- 
tal corrected blue magnitude, m Bcorr < 15.5 mag, and blue 
light isophotal angular diameter D25 > 1', while \b\ > 30°. 
The S 4 G data are probing down to stellar surface densities 
< 1 M pc~ 2 . The azimuthally averaged surface brightness 
profiles typically trace isophotes down to /i3.6^m > 27 (AB) 
mag arcsec' 2 . All of the images are processed using the S 4 G 
pipeline (see details in Sheth et al. 2010). The resulting pixel 
scale of S 4 G image is 0.75"/pixel and the typical full width 
half maximum of the Point Spread Function (PSF) is 1.7". 
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TABLE 1 

Catalog of parent sample of early type galaxies 



Object 


T-type 


Morphology 


D25 


B T 
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(4) 
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(6) 
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1.35 
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1.05 
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—16.69 
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JNCjC 0254 
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12.59 


—19. 1 1 


1624 
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1 i A 
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—18.73 


175 1 
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— 2.U 
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2.63 


12.38 


—20.43 


2372 
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—4.0 


SAB0~ 
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1 1 ii 
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iNLrC 0680 
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1 i on 
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mo 

lily 
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SAO 
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/l /IT 

4.4 7 
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1 1 1 c 
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1.48 
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JNCjC 1332 
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1524 
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NCjC 30/3 


—2.8 


SO 


1.20 


14.1 1 


-17.50 


1217 


NCjC 3384 


—2.7 




5.25 


10.89 


-19.82 


903 
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12.53 
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NGC 5338 


-2.0 
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14.05 


-16.61 
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NGC 548 1 


-3.9 




1.78 


13.28 


-19.26 


1883 


NGC 5574 


-2.8 




1.32 


13.29 


-18.79 


1657 


NGC 5576 


-4.8 




2.82 


11.79 


-20.14 


1482 


NGC 5846 


-4.7 


E+O/SAO" 


4.27 


11.09 


-21.27 


1824 


NGC 5854 


-1.1 




3.02 


12.65 


-19.70 


1750 


NGC 5864 


-1.7 




2.51 


12.68 


-19.76 


1850 


NGC 60 14 


-1.9 




1.74 


13.64 


-19.46 


2491 


NGC 6278 


-1.9 




1.74 


13.74 


-19.73 


2776 


NGC 7465 


-1.9 




1.07 


13.36 


-19.33 
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PGC 014037 


-1.0 




1.20 


15.12 


-17.74 


2546 


UGC 09348 


-1.6 




1.78 


14.72 


-17.45 
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NOTE. — Information on galaxies checked for tidal debris. All data were taken from Hyperleda, 
unless noted otherwise. Col. (1) : Object. Col. (2) : Numerical morphological type. Col. (3) : Galaxy 
morphology at 3.6/mj (Buta et al. 2010). Col. (4) : Projected major axis of a galaxy at the isophotal 
level of B = 25 mag/arcsec 2 . Col. (5) : Apparent total B magnitude. Col. (6) : Absolute B magnitude. 
Col. (7) : Heliocentric radial velocity obtained from radio observations. 



a : Heliocentric velocity from Smith et al. (2000). This galaxy does not have a value for v md in the 
Hyperleda database at the time of writing, hence it is not included in the S 4 G sample. But as the nearby 
galaxy, NGC 2634A is in the S 4 G sample, NGC 2634 was imaged in the same frame. 
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Using the zero magnitude flux at 3.6 fim, 280.9 Jy, as indi- 
cated in the IRAC Instrumental Handbook, 3.6 fim AB mag- 
nitude can be converted from 3.6 fim Vega magnitude using 
the following equation 



M 3 . 6tAB = Af 3 .6,v ega + 2.78 



(1) 



To identify tidal features in ETGs, we first chose ETGs 
from S 4 G using their numerical Hubble type T (de Vau- 
couleurs et al. 1992), which was obtained from Hyperleda 
(Paturel et al. 2003). We selected ETGs with T less than 
or equal to 0, thus our sample includes both Ellipticals and 
SOs. Note that although S 4 G surveys volume-, magnitude- 
, and size-limited samples, there is a caveat especially on the 
completeness of ETG samples. Since S 4 G chose galaxies with 
radial velocity Vhi < 3,000 km s -1 for sample selection, S 4 G 
misses gas-poor faint and small ETGs. We draw and analyze 
an initial set of 65 ETGs that were available at the moment 
of this pilot study, of which 34 are from the cryogenic mis- 
sion archival data and 3 1 are newly obtained during the post- 
cryogenic phase of Spitzer. The basic properties of the 65 
ETGs are listed in Table 1. Our sample consists of 16 ellipti- 
cals (T < -3.5; criteria adopted from Paturel et al. 2003 and 
Emsellem et al. 201 1) and 49 SOs (-3.5 > T < 0). There will 
be 180 ETGs in the complete sample of 2,300 S 4 G galaxies. 
ETGs with tidal debris presented in this study are located at 
distance between 20 and 43 Mpc, so the typical PSF (1.7" ) 
corresponds to a linear scale of 170 to 350 pc. 

2.2. Data Analysis 

The principal difficulty in identifying tidal debris is the very 
high contrast between the bright component of a galaxy and 
the faint, fine-scale tidal debris structures. To mitigate this 
contrast, we used the following methods: unsharp masking 
(Malin 1977), the building of a structure map (Pogge & Mar- 
tini 2002) and color map, and the two-dimensional profile fit- 
ting of the light distribution of the galaxy. Applying unsharp 
masking is a fast and easy way to remove the light associ- 
ated with the galaxy. It enhances the features of galaxies that 
exhibit sudden changes in their light distributions, which can 
be caused by spiral arms, shells and ripples. To generate an 
unsharp masked image, we first convolved the image with a 
circular Gaussian kernel of a = 15 pixels (corresponding to 
11.25"), a size that was chosen iteratively to better reveal 
the tidal debris in the unsharp masked images. We then di- 
vided the original image by the convolved image, which al- 
lowed us to highlight the sharply varying features from the 
smoother background. This technique was sufficient for re- 
vealing shells, arc-like features (see Figure 6), and highly- 
inclined disk components. 

Making structure maps (Pogge & Martini 2002) is another 
method to enhance unresolved or marginally resolved features 
on the scale of the PSF. This technique involves removing the 
large scale smoothed light of the galaxy using the image PSF 
as a kernel. The structure map therefore has the extra benefit 
of reducing the time spent in choosing the best kernel to re- 
veal the tidal features of the galaxy. To derive the structure 
map, we first divide the original image by the PSF-convolved 
image, then convolve this resulting image with the transpose 
of the PSF image. This operation can be expressed as follows, 



(2) 



I®P 



transpose of PSF image, and S is the image of the structure 
map (Pogge & Martini 2002). Smooth changes in the light 
distribution, such as faint tidal tails, are still not easily re- 
vealed with an unsharp masked image or structure map. 

Modeling the light distribution of the galaxy using its 
known components such as a bulge and/or a disk component 
is another common method, used for detecting tidal features. 
The Ellipse task in iraf can be used for modeling the light 
distribution of the galaxy. This task creates series of concen- 
tric elliptical isophotes that are fit to the galaxy image. The 
resulting galaxy model is subtracted from the original image 
to reveal the tidal features. This technique is, however, not op- 
timal to identify shells because concentric arc-like shell fea- 
tures can become incorporated in the ELLIPSE model. An- 
other problem with the ellipse-fitting method is that it can pro- 
duce spurious residuals when there are rapid changes in the 
ellipticity and position angle (PA). In such cases, the model 
subtraction leaves an imprint, which resembles shell-like de- 
bris (Janowiecki et al. 2010). We also found in our analysis 
that the residual images of the ELLIPSE model fit can show 
"butterfly" (quadrupole) patterns in the central region (also in 
Janowiecki et al. 2010). This occurs because the ellipse fit- 
ting cannot account for high order deviations from an ellipses 
(e.g., boxy bulges), which may be present near galactic cen- 
ters. 

A better technique for identifying tidal features is to use 
a multi-component, two-dimensional fitting algorithm like 
GALFIT (version 3.0, Peng et al. 2010; Peng et al. 2002). 
In order to find tidal debris, we examine the residual im- 
age, which were obtained by subtracting the best fit model 
from the image. To find the best fit model for each galaxy, 
we run GALFIT with four different functions; two different 
one-component fits that fit the bulge only and two different 
two-component fits that fit the bulge and the disk. These four 
functions are de Vaucouleurs, Sersic, Sersic+exponential and 
Sersic+Sersic profile. These fits will be discussed in more de- 
tail in Section 5.1. Although it has been shown that most of 
the ETGs do not follow a R l/4 law (Caon et al. 1993; Gra- 
ham et al. 1996), i.e., the Sersic index of the bulge is rarely 4, 
we modeled the ETGs with a n=4, de Vaucouleurs' profile to 
compare our results with previous studies on the KR and FP. 
Since some elliptical galaxies are mistaken SOs, we include a 
disk component for modeling those galaxies to increase our 
chances of diagnosing such cases. If we choose an unrealistic 
Sersic index to model galaxies, spurious features can also be 
produced in the residual image, therefore careful decomposi- 
tion is needed. 

The Sersic profile (Sersic 1963) is defined as, 



I(r) = 4ff exp < -b n 



Ceff) 



1/n 



1 



(3) 



where r e g is the effective radius, 7 e ff is the intensity at r = r e ff, 
n is the Sersic index, and b„ is chosen to satisfy 



/ I(r)2nrdr = 2 / l(r)2-nrdr. 
Jo Jo 



(4) 



where / is the original image, P is the PSF image, P' is the 



When n = 4, it is called the de Vaucouleurs profile, and if n = 1 , 
it describes the exponential disk profile. The scale length h r 
of the disk is given as, 

r eff =1.678ft r . (5) 
GALFIT gives the effective semi major axis in the output, 
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TABLE 2 

Basic Information of the Galaxies showing tidal features. 



Object 


Morphology 


Morphology 


Kinematic class 




Distance B To , a l Velocity dispersion 




(RC3) 


(Buta) 


(ATLAS 10 ) 


[deg] 


[Mpc] 


[Mag] 




[km s 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 




(8) 




NGC 0474 


SA(s)0 D 


(R)SAB0/a (shells) pec 


Fast Rotator 


31.3 


27.7 


12.38 


163 


± 


5<a) 


NGC 0680 


E + pec: 




Fast Rotator 


22.7 


34.8 


12.90 


214 


± 


6 <a) 


NGC 1222 


SO" pec: 




Slow Rotator 


72.7 


31.2 


13.15 


119 


± 


13(6) 


NGC 1482 


SA0 + pec sp 


Sa: sp 






25.0 


13.10 








NGC 2634 


El: 


SA(nl)0" (shells) 






31.7 


12.91 


191 


± 


6 (a) 


NGC 3032 


SAB(r)0° 




Fast Rotator 


0.8 


25.0 


13.07 


82 


± 


16*' 


NGC 4106 


SB(s)0+ 








34.3 


12.28 


176 


± 


17 W 


NGC 4694 


SB0 pec 




Fast Rotator 


2.0 


20.4 


12.29 


61 


± 


5(6) 


NGC 5018 


E3: 


SAB0~ (shells) pec 






43.3 


11.69 


211 


± 


17*' 


NGC 5574 


SB0- :? sp 




Fast Rotator 


4.8 


26.3 


13.29 


77 


± 


3 (/) 


NGC 5576 


E3 




Slow Rotator 


7.4 


23.9 


11.79 


164 


± 


13*' 



NOTE. — Col. (1): Object. Col. (2): Galaxy morphology in the optical blue band from RC3 catalog. Col. (3): Galaxy 
morphology at 3.6/^m (Buta et al. 2010). Col. (4): Kinematic class of galaxies from ATLAS^ D study (Emsellem et al. 201 1). 
Col. (5): Kinematic misalignment angle in degrees from ATLAS^ D study (Krajnovic et al. 2011), where sin \T/ =| sin(PA /7 / wf - 
PAi,>j) |. Col. (6): Luminosity distances which were corrected to the reference frame defined by the 3K microwave background 
radiation, from the NASA/IPAC Extragalactic database (NED). Col. (7): Total B magnitude, from Hyperleda. Col. (8): Central 
velocity dispersion, from (a): Simien & Prugniel (2000); (b): Table 4 of Wegner et al. (2003); (c): Dalle Ore et al. (1991); (d): 
Dressier et al. (1991) (f): Simien & Prugniel (2002) 



whereas other studies often use, the circularized effective ra- 
dius, defined as 

fb 



fleff=?"eff x \ -. 

V a 



(6) 



We also derive R ef f for comparison with other studies. Mean 
surface brightnesses within r = r e g are calculated with 



< Meff >= m buige + 2.5log[irR eS ], 



(7) 



where mbui ge is the bulge magnitude within the effective ra- 
dius, (b/a) is the axis ratio of the bulge component, and r e g is 
effective radius of the bulge. 

We constrained GALFIT to use the same central posi- 
tion of the galaxy for both the bulge and disk components. 
GALFIT (version 3.0 or higher) allows us to fit asymmetric 
features with Fourier modes and coordinate rotations. Un- 
like Ellipse in iraf, GALFIT does not build models using 
isophotal fits, but instead fits the given image with a suitable 
model all at once. Thus, it does not allow radial variations of 
PA and ellipticity but returns one global value of PA and el- 
lipticity for each component of the model. We estimated the 
sky background by selecting 10 different regions outside the 
galaxy, avoiding any obvious faint tidal features in the back- 
ground. We fixed the sky value while running GALFIT. All 
of the components were convolved with the PSF for IRAC 
3.6 /im-band, which was created from the sub-sampled IRAC 
Point Response Function (PRF) provided by the Spitzer Sci- 
ence Center 19 . Because images of the S 4 G are mosaicked by 
drizzling, and have a pixel scale of 0.75", we re-sampled the 
PRF 61 times (=0.75/1.22x100) to match the pixel scale of 
the S 4 G data. We also rotated each PSF to the average angle 
at which a galaxy was observed. For the Galactic extinction 
correction, we use the result of Schlegel et al. (1998). Be- 
cause the difference of the Galactic extinction between 3.6 
and 3.8 fim is not significant, we adopted the value of L' (3.8 
fj,m) for each object, taken from the NASA/IPAC Extragalac- 
tic Database (NED) 20 . 

19 http://ssc.spitzer.caltech.edu/irac/calibrationfiles/psfprf 

20 http://ned.ipac.caltech.edu 



To find the best fit model unaffected by bright foreground 
stars or faint background galaxies, we mask out those objects 
before running GALFIT. To generate mask images, we run 
Source Extractor (Bertin & Arnouts 1996) with a 3 sigma de- 
tection threshold and choose "check-image" to return "seg- 
mentation image" and make use of those as mask images. 
However, the "segmentation image" does not generate ade- 
quate masks, especially in the central parts of some galax- 
ies. So for sources lying within 0.5x£>25 from the center 
of galaxies, we create circular masks with the size propor- 
tional to the magnitude of the source. Then, we perform a 
preliminary GALFIT run to identify tidal features. Based 
on this first analysis, we improve the mask images manually 
by adding shells or diffuse tidal features to the masks, which 
do not originally include those features in the "segmentation 
image". Then, we run GALFIT again on this star-free, back- 
ground galaxy-free, and tidal feature-free image to find the 
best fit. In some cases, around the edge of masks, light from 
stars appears to be leaking out as can be seen in Figure 1 to 5. 
However, these images are displayed with a histogram equal- 
ized stretch that emphasizes faint structures and the leaks do 
not affect our results. 

3. GALAXIES EXHIBITING TIDAL FEATURES 

Of the 65 ETGs in our sample, 1 1 galaxies show evidence 
of tidal interactions, such as shells, ripples or tails in the 
GALFIT residual image. We also check the unsharp masked 
image, structure map. Especially for the diffuse features, we 
visually cross-check those features in the 3.6 fim image using 
the histogram stretch and in the GALFIT residual images. 
The basic properties of each galaxy are summarized in Ta- 
ble 2. We describe properties of galaxies showing shells and 
tidal debris in Section 3.1 and 3.2, respectively. All of 11 
galaxies were previously classified as galaxies with shells or 
tidal features, therefore we briefly summarize published stud- 
ies on them. Then we present the fraction of ETGs with shells 
and tidal disturbances and we compare them with literatures 
in Section 3.3. We also discuss the ring-like feature in NGC 
2974, which may not be tidal debris. 

3.1. Shell galaxies 
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FIG. 1. — Residual images of NGC 474 (left) and NGC 2634 (right). The model fit of two Sersic components, one for the bulge and one for the disk, is 
subtracted from each image. Left: Shells and tidal features of NGC 474. The"comma-like" feature is marked with 'a' atr~3.5'. Right: Shells are marked with 
'a' (r~45"), 'b','c' (r~75"), 'd' (r~100"), and a possible lens is marked with 'e' (r~25"). North is up and East is to the left for all of the images shown in this 



• NGC 474: This galaxy is also known as Arp 227 (Arp 
1966), and is interacting with NGC 470, which is 5.5' (~ 50 
kpc) to the West. An HI gas bridge connects NGC 474 and 
NGC 470 (Rampazzo et al. 2006). NGC 474 is famous for its 
shell structures (Malin & Carter 1983; Turnbull et al. 1999) 
and the misalignment between the kinematic and photometric 
major axis (Hau et al. 1996). About 3.5' to the South-west 
of the galaxy center, there is a "comma-like" structure (Ram- 
pazzo et al. 2006), which is denoted with a letter 'a' in Figure 
1. 

In the S 4 G data, shells in NGC 474 are clearly visible in the 
unsharp masked image and in the residual image (Figure 2) 
from the GALFIT model. There is also a possible bar feature 
(PA 33° from North to East), which shows up in the resid- 
ual image after a Sersic plus the exponential disk model has 
been subtracted (Figure 1). This is consistent with the findings 
of Laurikainen et al. (2010), who also noted that this galaxy 
has a weak oval bar surrounded by a lens. A lens is defined 
as a oval-shaped structure in the disk and its surface bright- 
ness is slowly decreasing with a sharp outer edge (Kormendy 
1979; Athanassoula 1983; Bosma 1983; Kormendy & Kenni- 
cutt 2004). We also find that shells and tidal features in the 
residual images at 3.6 /J,m are similar to those in R-band im- 
ages obtained by Turnbull et al. (1999). The PA and ellipticity 
of NGC 474 are in agreement with what Turnbull et al. (1999) 
have found. The 3.6-4.5 \im color map is shown in Figure 2. 
Morphological features of NGC 474 seen in the 3.6-4.5 /im 
color map are similar to those in the 3.6 [im image, and this is 
same as for other galaxies showing tidal features in this study. 

• NGC 2634: This galaxy presents multiple shells at r <~ 45, 
75 and 100", which correspond to projected distances of r ~ 7 
kpc (marked with 'a' in Figure 1), 1 1 kpc ('b' and 'c'), and 15 
kpc ('d')- In the residual images (Figure 3), a relatively bright 
ring-like feature is present at r~25" as marked with 'e'. It 
is difficult to decide whether this is a shell or not because it 
does not show a sharp edge on the surface brightness profile 
and it encircles the central part of the galaxy. Therefore, it is 
more probable that the ring-like feature is a lens, as Buta et al. 
(2010) classified this galaxy has a nuclear lens. There are only 



a few previous studies that have noted these shells. While 
the galaxy was originally classified as "El:" in the RC3, it 
was recently classified by Buta et al. (2010) as a SA(nl)Q~ 
(shells) from the S 4 G data. In a conference abstract, Statler 
et al. (2001) report that they found shells in NGC 2634 using 
optical image. 

• NGC 3032: The galaxy has shell at r-30" and it is 
marked with arrows in Figure 7. In the residual image, we 
see some features in the central region. We check the F606W 
HST and Sloan Digital Sky Survey (SDSS) image and found 
that NGC 3032 has central dusty spiral features within 11". 
Schweizer & Seitzer (1988) also showed that this galaxies 
have 3 interleaving ripples in Figure 2 of their paper. Lau- 
rikainen et al. (2010) noted that this galaxy has two nearly 
circular lenses extending out to r^7 and 16", respectively. 

• NGC 5018: NGC 5018 shows multiple shells (Malin & 
Carter 1983) and broad fans of stellar light (Figure 4). One of 
the inner shells is in the shape of an hourglass, which is also 
evident in the V-band image in Figure 2 of Fort et al. (1986). 
Dust lanes and patches in the galaxy have been noted (Fort 
et al. 1986; Sandage & Bedke 1994) and with a UV F336W 
HST image, Rampazzo et al. (2007) find dust lanes also in 
the nuclear region ( 1 0" x 1 0") of the galaxy. A filament of HI 
gas stretches across the stellar body of NGC 5018 and con- 
nects two neighbors (Serra & Oosterloo 2010). A faint dis- 
persed tidal tail lies toward the North-west. The shells are not 
detected in a near-ultraviolet image taken from the GALEX 
archive (Rampazzo et al. 2007), indicating that there is no 
strong recent star formation in shells. The unsharp masked 
image, structure map, and residual images (Figure 4) suggest 
that there is an inner disk inclined at a PA of 97° in the central 
region of the galaxy, extending out to 25" from the galactic 
center. 

3.2. Galaxies showing tidal features 

• NGC 680: This galaxy is classified as a "E + pec:" (RC3), 
and it belongs to the NGC 691 group, (van Moorsel 1988; 
Crook et al. 2007) with NGC 678 (at a projected distance of 
5.4', 60 kpc ), NGC 694 (16.6', 190 kpc), NGC 691 (17.9', 
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FIG. 2. — Image analysis of NGC 474 and NGC 680. TopQeft to right): 3.6fim image; unsharp masked image; structure map; and color image of 3.6-4.5 fim 
in histogram equalized stretch. Bottom (left to right): GALFIT residual images for a bulge (n=free) and an exponential disk (n=l) fit; de Vauouleurs law (n=4) 
fit; bulge-only (n=free) fit, and bulge (n=free) with disk (n=free) fit, respectively. The white circle on the bottom left of the unsharp masked image indicates the 
size of the Gaussian kernel, which was used to generate the unsharp masked image. The white circle in the residual image of the 'Ser+Exp' model presents the 
adopted inner boundary for estimating the fractional luminosity of tidal debris, which is also listed in Table 4. 



200 kpc, also known as NGC 674), and IC 167 (19.1', 220 
kpc). In our analysis of the S 4 G data, we clearly see tidal dis- 
turbances, shells, and an asymmetric envelope in the residual 
maps from GALFIT shown in Figure 2. As an exercise, we 
have also run the ELLIPSE task in IRAF with a varying cen- 
ter to check if we can confirm the presence of the asymmetry 
in NGC 680. With this fit we can reproduce the asymmet- 
ric envelope. However, the model image produced from the 
ELLIPSE fit shows artificial structures caused by the pile up 
of ellipses. Using this model to produce residual image simi- 
larly as with our GALFIT models, we still find the same tidal 
features, although their shape might vary. The dark feature at 
the South-East side of NGC 680 in Figure 2 becomes smaller 
in size - and part of it becomes bright - but the bright feature 
at the North- West, which we describe as a shell, still remains. 
NGC 680 has also been described to have possible patchy dust 
(Ebneter et al. 1988). Due et al. (201 1) showed that NGC 680 
exhibits two diffuse plumes in the optical image and these are 
connected to an extended HI tail. 

• NGC 1222: There is a tidal disturbance toward the North 
as can be seen in Figure 3. Also known as Markarian 604, 
NGC 1222 has a starburst nucleus, which has a compact radio 
and infrared source (Beck et al. 2007). Petrosian & Burenkov 



(1993) claimed that this galaxy has appearance of a triple sys- 
tem in Figure 1 of their paper, and argued that two other galax- 
ies are dwarf galaxies with high surface brightness interacting 
with the main galaxy. They also mentioned that there is a con- 
tinuous twist in the isophotes. To check if we can confirm this 
isophotal twist, we have run the ELLIPSE task in IRAF, and we 
also find the twist in the isophotes of the galaxy. Also NGC 
1222 have a large kinematic misalignment between the pho- 
tometric and kinematic major axes (^ = 12.1°, where sin W =| 
sin(PAp/, of - PA&„) |). Despite the complex morphology, NGC 
1222 is classified as a slow rotator (Emsellem et al. 201 1). 

• NGC 1482: Toward the Northeast, the galaxy has a tidal 
feature, and it is marked with an arrow in Figure 7. Hameed 
& Devereux (1999) note that this galaxy has "filaments and/or 
chimneys of ionized gas extending perpendicular to the disk" 
(see Figure 5 of their paper). Large scale hourglass-shaped 
outflows has been detected with excitation maps of [N II] 
A6583 and Ha also by Veilleux & Rupke (2002). They 
claimed this galactic wind is extending along the minor axis 
of the galaxy at least 1.5 kpc above and below the galactic 
plane in both directions with a velocity of ^250 km s" 1 . 

• NGC 4105/4106 pair: These galaxies are interacting and 
NGC 4106 (left one in Figure 7) shows tidal disturbances. 




Fig. 3.— 



Image analysis of NGC 1222, NGC 2634 and NGC 4694, same as figure 2. 
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Fig. 4.— 



Image analysis of NGC 5018 and NGC 5576, same as Figure 2. 
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NGC 5576 / NGC 5574 



2 arcmin 





FIG. 5. — The upper left panel shows the 3.6 \xm S 4 G image of NGC 5576 (upper left galaxy) and NGC 5574 (lower right galaxy). The upper right panel 
shows the GALFIT residual image of the model fit with Sersic+exponential, displayed in histogram equalized stretch. A long tidal tail extends out to 7' and 
tidal disturbances of NGC 5576 are also present. The companion galaxy, NGC 5574, and stars were masked out to enhance faint structures. The lower left 
panel presents the regions that were chosen to estimate the sky brightness value in the NGC 5576 field at 3.6 /im, displayed in a histogram equalized stretch with 
foreground stars masked out. The lower right panel shows the regions that were adopted to compare the brightness of the tidal tail with that of the sky. 



One of the tidal tails extends toward the elliptical galaxy, 
NGC 4105 (right one in Figure 7). These galaxies are sep- 
arated by a projected distance of 7.7 kpc and the difference in 
radial velocities is ^260 km s" 1 . NGC 4106 shows a disky 
central structure and NGC 4105 has dusty features in the cen- 
tral region and disky outer structure (Griitzbauch et al. 2007). 
Reduzzi & Rampazzo (1996) obtained a nearly constant color 
profile of (B-V) ~ 0.9 for both NGC 4105 and NGC 4106, 
while Longhetti et al. (1999) claimed that these galaxies expe- 
rienced recent star formation episodes from their study of line 
indices. While NGC 4106 is perturbed, NGC 4105 does not 
show clear perturbation (Griitzbauch et al. 2007). Thus, we 
exclude the elliptical galaxy, NGC 4105, from our selection 
of galaxies with tidal disturbances. 

• NGC 4694: This galaxy is described as "amorphous" 
(Sandage & Tammann 1981), and classified as an SBO(pec) 
in RC3 (de Vaucouleurs et al. 1991). However, as can be seen 
in Figure 3, it is not easy to tell whether NGC 4694 exhibits 



any bar feature in the 3.6 /im image. Using multi component 
galaxy decomposition, Laurikainen et al. (2010) note that it 
has a bright oval-like component and a peculiar nuclear bar or 
an inner disc. Due et al. (2007) presented a VLA HI map of 
NGC 4694 and its nearby dwarf galaxy, VCC 2062, in their 
Figure 1 . They show that the two galaxies are connected by 
an HI bridge. The faint spiral arm-like tidal feature of NGC 
4694 in the bottom panel of Figure 3 likely follows Hi-arm. 
Therefore, this tidal feature may originate from an interaction 
with VCC 2062. 

• NGC 5576: NGC 5576 is reported to have boxy isophotes 
(Ebneter et al. 1988), and it shows peculiar tidal features, 
tidal disturbances to the North and a very long tidal stream 
to the South. The tidal disturbance in this galaxy can be seen 
in the right panel of Figure 5. At an angular separation of 
about 2.76' (~ 16 kpc) to the South-west from NGC 5576 
(V rarf = 1482 km s" 1 ) there is a barred SO type galaxy, NGC 
5574 (N ra d = 1657 km s" 1 ). There is a long tidal tail extend- 
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FIG. 6.— Left: NGC 2974 at 3.6 pm Right: Unsharp masked image of NGC 2974. 




FIG. 7. — Image analysis of NGC 1482, NGC 3032 and NGC 4106, which were identified to have tidal features. Note that these galaxies were not carried on 
for the further analysis, because of the muxbleed (marked with 'a' in the upper panel) and column pull-down (marked with 'b') due to saturation (NGC 1482 and 
NGC 3032) or close companion (NGC 4105 and NGC 4106). (Left to right): 3.6 fan image in log-scale stretch; 3.6 jim image in histogram equalized stretch; 
structure map; GALFIT residual images from a bulge (n=free) and an exponential disk (n=l) fit. Tidal features are marked with arrows. 
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ing out to 60 kpc, which suggests that NGC 5576 and NGC 
5574 are an interacting pair. To check how bright this tidal 
tail is compared to the sky, we chose square regions outside 
the galaxy, as shown in the bottom left panel of Figure 5. We 
then estimated the surface brightness of the tidal tail in these 
selected boxes. The signal to noise ratio of the tidal tail re- 
gion is low (S/N ~1.5), which indicates that the tidal tail is 
only marginally brighter than the sky. Tal et al. (2009) also 
see this faint tidal tail in their deep V-band image out to 75 
kpc. 

Despite the faintness of the tidal tail, NGC 5576 and NGC 
5574 are likely interacting. The nucleus of NGC 5574 has 
an iron-overabundant stellar system (Sil'chenko et al. 2002), 
which is a rare phenomenon for ellipticals, because most ellip- 
ticals are magnesium-overabundant. Sil'chenko et al. (2002) 
mention that this iron overabundance can be observed in some 
irregular galaxies and is usually interpreted as a result of mul- 
tiple starbursts. Interestingly, the nucleus of the companion 
galaxy, NGC 5576, also shows this rare iron overabundance. 
The luminosity-weighted stellar ages of both galaxies are 3 
Gyr (Sil'chenko 1997), which implies that these two galaxies 
have experienced similar recent episodes of star formation. 
NGC 5576 is classified as a slow rotator (Emsellem et al. 
2011) without any specific kinematic feature, such as kine- 
matically distinct core or kinematic twist in the central part 
of the galaxy. This result implies that this galaxy experienced 
a significant merging event without being able to rebuild a 
fast rotating stellar component through dissipation (Emsellem 
et al. 2011). Although the long tidal tail seems to be related 
to the interaction of NGC 5576 and NGC 5574, and even if 
NGC 5576 shows tidal disturbances in the 3.6 fim imaging 
data alone, we can not be sure from which galaxy the tidal tail 
originated. As catalogued in the Isolated triplets of galaxies 
(Karachentseva et al. 1979, KTG 55), NGC 5576 has another 
companion galaxy, NGC 5577, located to the North, which 
has a redshift similar to NGC 5576. Since we only cover 
the southern part of this galaxy pair at 3.6/im, we checked 
the northern part in the accompanying 4.5 [im data, but could 
not identify a tidal bridge between NGC 5576 & NGC 5577. 
Therefore, the tidal disturbance of NGC 5576 might be related 
to the interaction with NGC 5574, although NGC 5574 has a 
larger recession velocity or it may be that NGC 5574 has a 
high peculiar velocity due to its interaction with NGC 5576. 

• NGC 2974: This object is classified as an elliptical galaxy, 
E4 (de Vaucouleurs et al. 1991), isolated elliptical (Sandage 
& Bedke 1994), and an SA(r)0/a galaxy by Buta et al. (2010) 
using the S 4 G data. Tal et al. (2009) identified multiple shells 
in this system. NGC 2974, however, shows an HI velocity 
field consistent with a rotating disk and the rotation axis of the 
galaxy is well aligned with the optical minor axis of the galaxy 
(Kim et al. 1988). Additionally, according to Emsellem et al. 
(2003), this galaxy has gaseous spiral structure and a possi- 
ble inner bar within a radius of 200 pc. From the near- and 
far-ultraviolet study (Jeong et al. 2007), NGC 2974 hosts an 
outer ring at a radius of 60" with a young stellar population 
(< 1 Gyr), and we confirm the arc-like ring features at a radius 
of 55" (corresponds to 7 kpc) in the unsharp masked image of 
3.6 /im (Figure 6). A bright star near NGC 2974 makes it dif- 
ficult to determine whether this feature is a ring or a shell, but 
the arc-like feature encompasses almost 180°, which can be 
seen in Figure 1 of Jeong et al. (2007). It also hosts a regular 
HI ring aligned with the stellar component of galaxy (Serra & 
Oosterloo 2010). Therefore, this arc-like feature may not be 



a shell as previously claimed (Tal et al. 2009). For these rea- 
sons, we excluded NGC 2974 from our sample of ETGs with 
tidal features. 

3.3. Fraction of tidally disturbed ETGs 

We have identified 11 tidally disturbed galaxies from 65 
ETGs. Note that our sample of this pilot study does not in- 
clude the whole statistically complete sample of ETGs in S 4 G. 
We found fewer ETGs with tidal features than van Dokkum 
(2005) and Tal et al. (2009), who found fine structures in 70% 
of their ETGs. Tidal features can be relatively faint, as shown 
by Janowiecki et al. (2010), the peak surface brightness of 
tidal features range from fi v = 26 - 29 mag arcsec" 2 . There- 
fore, the ability to identify tidal features depends strongly on 
the data, especially on the depth of the image. We summa- 
rize the fraction of tidally disturbed galaxies from the liter- 
ature in Table 3 along with the depth of images analyzed in 
each study. As studies in the literature were performed in dif- 
ferent bands, we convert the magnitude of depth of images 
to B-band AB magnitude system in Table 3 to directly com- 
pare the fraction of tidally disturbed ETGs. As we can see 
in Table 3, the deeper studies can find more tidally disturbed 
ETGs. Depending on the project, the Spitzer archival images 
have different depths and we could find tidal features down 
to /i3.6,AB < 26.5 mag arcsec" 2 , while in the newly obtained 
"Warm" data from the S 4 G project, we found tidal features 
down to H3.6AB < 25.2 mag arcsec" 2 (<~ ^b,ab< 26.9 mag 
arcsec" 2 ). 

S 4 G maps each galaxy to 1.5x^25. However, for one of 
our sample galaxies, NGC 5846, which is identified to have 
shells by Tal et al. (2009), we could not identify those shells 
because they lie outside of our field. As discussed in Sec- 
tion 3.2, we exclude NGC 2974 when we count ETGs show- 
ing shells or tidal disturbances. IC 51 (Malin & Carter 1983) 
and NGC 3608 (Forbes & Thomson 1992) are found to have 
shells from the literature 21 . But, we did not find clear evi- 
dence for tidal features. With magnitude and field of view 
limit, we find 17±3% (11/65) of ETGs are tidally disturbed 
and 6±2% (4/65) have shell. If we restrict our sample to ellip- 
ticals using the morphology information of RC3, we find that 
13±1% (2/19) of ellipticals have shells to our limiting surface 
brightness and field of view. We make use of Poisson statis- 
tics to present an uncertainty estimate of our fraction, using 
a = [f(l - f/N)] x l 2 where N is the sample size and f is the 
number of galaxies showing tidal debris or shells following 
the method of Laine et al. (2002). 

4. LUMINOSITY OF TIDAL FEATURES 

To estimate the luminosity of the tidal features, we first 
mask out all the foreground stars and background galaxies, 
and after checking the residual image of each galaxy, we also 
mask out known features such as lens components or an in- 
ner bar, mentioned in literature. The masked out inner and 
outer boundaries are listed in Table 4. Because now we focus 
on the tidal features, this time we do not mask out tidal fea- 
tures but leave them in the GALFIT residual image. Then, 
we estimate the ratio of the total counts from the tidal debris 
in the residual image to the total counts in the model image. 

21 NGC 4494 is in the sample of shell galaxies in Table 4 of Forbes et al. 
(1994), but afterwards there was no more description on the shell in this 
galaxy (Forbes et al. 1995a, Foster et al. 2011). So we exclude NGC 4494 
from the sample that have been reported to have a shell in the optical. 



ETGs with Tidal Debris and their Scaling Relations in S 4 G 



13 



TABLE 3 

Fraction of ETGs showing tidal features. 



Study 


Fraction 


Object 


Band 


Depth 


Equivalent depth in B 




[%] 






[mag] 


[AB mag] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Malin & Carter (1983) 


IT (shell) 


Ellipticals 


B 


26.5 (Vega) 


26.3 


van Dokkum (2005) 


70 


ETGs 


BVR b 


29.5 b (AB) -> 28.0 C 


29.5 b ->■ 28.0 C 








BVI b 


29.0 b (AB) -> 28.0 C 


29.0 b ->■ 28.0 C 


Tal et al. (2009) 


73 


ETGs 


V 


29.0 (Vega) -» 27.7 C 


30.0 -> 28.7 C 


Krajnovic et al. (2011) 


8 


ETGs 


r 


26.0 (AB) 


27.3 


This study 


17±3 


ETGs 


3.6 fim 


27.0 (AB) (25.2f 


28.7 -> (26.9) e 




13 a ±l (shell) 


Ellipticals 










6 f ±2 (shell) 


ETGs 









NOTE. — Col. (1): Author of the study. Col. (2): Fraction of tidally disturbed galaxies. Col. (3): Morphology of parent sample. Col. 
(4): Filter. Col. (5): Typical 1 a surface brightness limit presented in the literature. Col. (6): Typical 1 a surface brightness limit converted 
to B-band AB magnitude. We made use of the Table 13 of Eisenhardt et al. (2007) to convert magnitudes, and we took colors (B - V) = 1.0, 
(B-R) = 1.5 at H=13 mag. To convert r to R, we used the equation listed in SDSS web page. d We took the value (g-r) = 0.7 for ETGs and 
we used (R-r) = -0.2 for conversion. We adopt (3.6 fim - B) ~ 1.7 as used in Sheth et al. (2010). 
a Fraction of shell galaxies in ellipticals galaxies. 
b Co-added. 

c Depth that van Dokkum (2005) and Tal et al. (2009) can confidently identify features, mentioned in their paper. 
^ http : // www.sdss.org/dr5 /algorithms / sdssUBVRITransform.html#Lupton2QQ5 

e The surface brightness of the tidal features in the shallowest image of our sample. With deeper images from the Spitzer Heritage Archive, 
we could identify tidal features down to /i3.6,AB = 26.5 (~ fiBAB = 28.2 ), see the text for detail. 



Fraction of shell galaxies in ETGs, including elliptical and SO galaxies. 

We calculated the luminosity fraction of the tidal features as 
a percentage, {F residua i / F mode i) x 100 and list these in Table 4. 
The luminosity fractions of tidal features are estimated both 
in a circular region and in an elliptical region that is chosen to 
have the same PA and ellipticity as the bulge component of the 
model has. The ratios using the circular regions are labelled 
as Ratiodrc, and those using the elliptical region are labelled 
as Ratio e ii. In the cases of the Sersic+exponential fit and the 
Sersic+Sersic fit, the fractional luminosity is estimated using 
the total model, which is the sum of both the bulge and the 
disk. For simplicity, we present these ratios in the first row of 
S+E and S+S in Table 4. 

Some of our sample galaxies (NGC 1482 and NGC 3032) 
are severely saturated and have the "muxbleed" and "column 
pull-down" artifacts which result from a bright source - see 
the IRAC Instrumental Handbook Section 7.2 for more de- 
tails. The "muxbleed" is a trail of bright pixels along a row in 
the readout direction that repeats with an increased signal ev- 
ery fourth column. This is denoted as 'a' in NGC 1482 of Fig- 
ure 7. The "column pull-down" artifact has a decreased sig- 
nal above and below the bright source along a column. This 
is denoted as 'b' in Figure 7. Note that we excluded galax- 
ies which were severely contaminated by these effects in our 
analysis. Close companions (NGC 4105/ 4106 pair) that lie 
within 1/4 of D25 were also excluded. 

If we include all the pixel values in the residual images to 
estimate total counts of tidal features, then we measure the lu- 
minosity fraction of the tidal feature to be 0.5 - 7% of the total 
model galaxy luminosity. But as described in Bennert et al. 
(2008), to avoid adding negative values from the regions that 
were oversubtracted by GALFIT in the residual image, now 
we only sum values of all pixels brighter than —la of the sky 
background to get total counts, and list luminosity fractions in 
Table 4. If we only sum values of pixels that are brighter than 
-3cr of the sky background, then we obtain the luminosity of 1 

- 8% of the total galaxy luminosity. For galaxies with shells, 
this result is consistent with those of Prieur (1988), Canalizo 
et al. (2007), and Bennert et al. (2008), who found that the 
shells constitute 5% of B-band, <~ 6% of F606W-band, and 5 

- 10% of F606 W-band luminosity of the galaxy, respectively. 



We also estimate the tidal parameter, which were used in van 
Dokkum (2005) and Tal et al. (2009). The tidal parameter is 
defined as, 

F galaxy = \^x,y ~ M X y | , (8) 

where I X y and M^ y are pixel values at the position x, y of the 
object and model, respectively. The tidal parameters of our 
sample galaxies range from 0.09 to 0.23. These values are 
comparable to what van Dokkum (2005) and Tal et al. (2009) 
found for galaxies that have tidal features. 

In Section 3.2, for the galaxies which show isophotal twist 
or asymmetric envelope, we have run the ELLIPSE task in 
IRAF to confirm the presence of those features. We also es- 
timate the luminosity of tidal features using the fit from EL- 
LIPSE. The luminosity of the tidal features does not change 
significantly in the case of NGC 680, but is higher in the case 
of NGC 1222 with ellipse fits. Because of the artificial fea- 
tures shown in the models derived from ELLIPSE fits, and for 
the sake of consistency, we will only consider here the results 
derived from the GALFIT models. 

By assuming the same stellar M/L in shells and galaxies, 
we estimate that 3% to 10% of the galaxy stellar mass is con- 
tained in the shells or tidal features. The amount of mass in 
shells can be used to infer the type of merger that led to its 
creation. Canalizo et al. (2007) discuss that with their own nu- 
merical simulations and those of Hernquist & Spergel (1992), 
if shell features were formed through major mergers, the stars 
in the resulting shells constitute one fourth or less of the to- 
tal mass of the companion galaxy. Following the argument of 
Canalizo et al. (2007), shells with 5% of the mass are likely 
to result from a merger with a 3:1 mass ratio, whereas shells 
with 10% of the mass are more likely to come from a merger 
where the galaxies have comparable masses. 

The luminosity of the long tidal stream of the NGC 
5576/5574 pair is estimated in the same manner. We find that 
the long tidal stream constitutes 4% of the galaxy model lu- 
minosity of NGC 5576. The estimated luminosity of the tidal 
features here might change slightly due to foreground stars or 
background galaxies that might not be fully masked. We in- 
crease the size of mask by 3 pixels for every foreground stars 
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and background galaxies, and we find that the amount that 
long tidal stream constitutes decreased to 3% of the galaxy 
model luminosity of NGC 5576. 

5. STRUCTURAL PROPERTIES 

5.1. Modeling light distribution of galaxies with different 
functions 

The structural properties such as sersic index, effective ra- 
dius, disk scale length, position angle, minor to major axis 
ratio, and bulge to total of ETGs having tidal features are 
presented in Table 4. As noted in Section 2.2, we perform 
GALFIT twice, the preliminary GALFIT run is on the image 
where foreground stars and background galaxies were masked 
out and the second GALFIT run is on the image where tidal 
features were also masked out. The results in Table 4 are ob- 
tained with the images where tidal features were also masked 
out. The differences between two different runs show that the 
Sersic indices, effective radii of bulges, and scale lengths of 
disks decrease by less than 10%, except for NGC 474, where 
the effective radius of the bulge component that was obtained 
from the Sersic only model fit decreases by ~ 20% from the 
first preliminary GALFIT run to the second GALFIT run. 

We exclude galaxies that have muxbleeds and close com- 
panions from our final sample as discussed in the previous 
section and present the results for seven ETGs. Because not 
all ETGs follow the R l l A law, we need to test several model 
profiles to see which one best describes the two-dimensional 
light distribution of these galaxies. Even though galaxies were 
once classified as ellipticals, some of those ellipticals might 
actually be misclassified SOs. By including a disk compo- 
nent in the model, we increase our chance of diagnosing such 
cases. Some of those examples are NGC 2634 and NGC 5018. 
They were first classified as ellipticals in the RC3 using the 
optical band, but recently they were reclassified as SOs (Buta 
et al. 2010) from 3.6 [im data. In fact, except for NGC 680 
and NGC 5576, all galaxies are classified as SOs, and there- 
fore, we also present disk scale lengths of these galaxies in 
Table 4. We see that parameters such as the Sersic index, 
the effective radius or the disk scale length change with cho- 
sen model. We compare the differences in the effective radii 
in Figure 8. The left panel shows differences in the effec- 
tive radii obtained from the Sersic model fit and the de Vau- 
couleurs R X I A model fit. For most galaxies, the effective radii 
obtained from the de Vaucouleurs R l l A model fit are smaller 
than those from the Sersic model fit and they cluster between 
2-3 kpc. The middle panel of Figure 8 compares the effective 
radii of the bulge components, estimated from a fit of Sersic 
(for bulge, n=free) + exponential disk model fit and Sersic (for 
bulge, n=free)+Sersic (for disk, n=free) model. For the Ser- 
sic (n=free)+Sersic (n=free) model, we adopt the component 
with a larger Sersic index as a bulge and the one with smaller 
Sersic index as a disk. For NGC 1222 and NGC 2634, the ef- 
fective radius of the larger Sersic index component (bulge) is 
larger than that of the smaller Sersic index component (disk). 
Therefore, for those two galaxies, this model profile fit may 
not be realistic, but for reference we keep those data in Figure 
8, 9, and 10. If we introduce a disk component to our model 
galaxies, the effective radius of the bulge component becomes 
smaller. We see this in the distribution of effective radius in 
the left and middle panel. Therefore, we must be careful not 
to blindly classify all galaxies as ETGs and treat them as a 
group, but we must model the lenticular galaxies including a 
disk component. 



If we let GALFIT model a disk component with a free Ser- 
sic index, we find that the Sersic index of the disk is never 
equal to 1, as would be expected for a conventional exponen- 
tial disk. This may be because disks of galaxies cannot be 
characterized by a single Sersic index (Freeman 1970). In spi- 
ral galaxies, the outer disks show truncations (down-bending) 
or anti-truncations (up-bending) profiles (e.g., Pohlen et al. 
2002, Erwin et al. 2005, Munoz-Mateos et al. 2009, Munoz- 
Mateos et al. 2012, in preparation). Even though the Sersic in- 
dex deviates from n=l, the effective radius of the exponential 
disk from the Sersic+exponential model fit and Sersic+Sersic 
model fit is quite similar, as seen in the right panel of Figure 
8. 

We tabulate kinematic class (Emsellem et al. 2011) and 
kinematic misalignment angle (Krajnovic et al. 2011) in Ta- 
ble 2. If galaxies show tidal features, there should be im- 
prints on the kinematics of those galaxies, and we check the 
result of our sample galaxies against those of the ATLAS 3D 
study. NGC 474, NGC 680 and NGC 1222 show a rather large 
kinematic misalignment (^ > 20°) between the photometric 
and kinematic major axes. Only 10% (25/260) of ATLAS 3D 
galaxies have > 20° and these are extreme cases. We find 
that one of our shell galaxies (NGC 474) and the galaxies with 
tidal disturbances (NGC 680 and NGC 4694) are better fitted 
with a disk component - these galaxies were classified as fast 
rotators (Emsellem et al. 201 1). This result is consistent with 
the finding of Cappellari et al. (2011b) that fast rotators are 
generally lenticular galaxies. One exception is NGC 5576. 
We find that this galaxy is better described with disk com- 
ponent, but it is classified as a slow rotator (Emsellem et al. 
201 1) and according to Krajnovic et al. (201 1), it is classified 
as a non-regular rotator, which does not show any specific fea- 
ture on the velocity map, such as kinematically distinct core 
or kinematic twist or counter rotating core. But Krajnovic 
et al. (2011) also note that it is possible that in lower noise 
velocity map, the galaxy would be classified as a regular rota- 
tor, which implies it might have a disk component. Thus our 
results from model fits are in agreement with the kinematic 
study from ATLAS 3D . 

We find that the total magnitudes that were obtained as- 
suming a specific model fit can vary among the different fits. 
For example, the total magnitude of NGC 474 varies up to 
0.3 magnitudes between n=4, de Vaucouleurs' fit and n=free 
Sersic fit. If a galaxy is better described by a n<4 fit, but 
modeled with n=4, then we obtain a brighter total magnitude. 
On the other hand, if a galaxy is better described by a n>4 fit, 
but modeled with n=4, then we obtain a fainter total magni- 
tude. For NGC 1222 and NGC 4694, Sersic indices of one- 
component Sersic fits are smaller than 4, and the total magni- 
tudes obtained with those Sersic models are larger than those 
from de Vaucouleurs' fits. For other galaxies that have n>4, 
the opposite holds. Studies of ETGs often estimate total mag- 
nitudes assuming their light profiles as de Vaucouleurs' law. 
But, as shown here, it is important to consider the differences 
due to the fitting model to obtain reliable estimates. 

5.2. Scaling relations 

Surface brightnesses and effective radii are determined by 
the choice of the adopted model profile, as evident in Fig- 
ure 8, 9 and 10. Therefore, we need to be careful when we 
compare our results with those from other studies, especially 
when we compare with galaxies that were modeled with dif- 
ferent types of fit. We present the KR, a relation between the 



ETGs with Tidal Debris and their Scaling Relations in S 4 G 



15 



TABLE 4 
Results of modeling galaxies 



Object Model 


n 


type 


''cff 


'eff 


h T 




PA 


b/a 


mag^.b,AB MagT,, 6AB 


Ratio c { rc 


Ratio e ii 


B/T 


x, 2 


[inner.outer] 






["] 


[kpc] 


["] 


[kpc] 


[°] 




[mag] 


[mag] 


(Circular) (Elliptical) 






(1) (2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 






U ■>> 


(It) 


(15) 


(16) 


NGC 0474 S 


6.65 


free 


31.78 


4.86 






11.3 


0.91 


10.88 


-21.65 


6.2 


6.3 




1.14 


[45.0,217.5] S+E(B) 


2.44 


free 


4.46 


0.68 






8.3 


0.90 


12.14 


-20.39 


11.0 


11.1 


0.39 


0.82 


S+E (D) 


1.00 


fix 


40.30 


6.16 


24.02 


3.67 


-15.2 


0.94 


11.67 


-20.85 










DeV 


4.00 


fix 


14.29 


2.19 






12.7 


0.93 


11.23 


-21.30 


19.1 


18.5 




1.90 


S+S (B) 


2.34 


free 


4.18 


0.64 






7.7 


0.90 


12.20 


-20.33 


10.4 


10.5 


0.37 


0.81 


S+S (D) 


1.12 


free 


39.81 


6.09 






-10.6 


0.94 


11.62 


-20.91 










NGC 0680 S 


5.76 


free 


15.93 


2.93 






-10.5 


0.78 


11.06 


-21.87 


5.5 


5.4 




2.44 


[12.8,127.5] S+E(B) 


3.38 


free 


4.91 


0.90 






-4.3 


0.78 


11.84 


-21.09 


6.7 


7.2 


0.55 2.06 


S+E (D) 


1.00 


fix 


32.17 


5.91 


19.17 


3.52 


-23.7 


0.75 


12.07 


-20.86 










DeV 


4.00 


fix 


10.55 


1.94 






-9.0 


0.78 


11.23 


-21.71 


12.2 


12.5 




3.75 


S+S (B) 


2.57 


free 


2.98 


0.55 






0.2 


0.77 


12.28 


-20.65 


5.7 


6.1 


0.36 


2.01 


S+S (D) 


1.63 


free 


26.47 


4.87 






-21.0 


0.76 


11.66 


-21.27 










NGC 1222 S 


3.25 


free 


5.82 


0.94 






-73.5 


0.68 


12.05 


-20.59 


14.6 


15.5 




3.84 


[24.6,120.0] S+E(B) 


2.01 


free 


3.73 


0.60 






-78.1 


0.59 


12.42 


-20.22 


7.3 


8.8 


0.66 2.58 


S+E (D) 


1.00 


fix 


28.01 


4.51 


16.69 


2.69 


-15.6 


0.61 


13.16 


-19.48 










DeV 


4.00 


fix 


6.70 


1.08 






-72.4 


0.69 


11.99 


-20.65 


11.5 


12.5 




4.06 


S+S (B) 


2.82 


free 


21.37 


3.44 






-19.4 


0.68 


12.63 


-20.01 


5.1 


6.2 


0.51 


2.41 


S+S (D) 


1.89 


free 


3.43 


0.55 






-80.6 


0.50 


12.69 


-19.95 










NGC 2634 S 


6.41 


free 


29.59 


4.40 






32.7 


0.91 


11.22 


-21.25 


3.9 


4.3 




1.23 


[12.8,142.5] S+E(B) 


4.06 


free 


9.25 


1.38 






29.1 


0.88 


11.96 


-20.50 


5.0 


6.1 


0.61 


1.17 


S+E (D) 


1.00 


fix 


47.61 


7.08 


28.38 


4.22 


74.7 


0.95 


12.47 


-20.00 










DeV 


4.00 


fix 


15.23 


2.26 






33.1 


0.91 


11.49 


-20.97 


13.7 


14.5 




1.92 


S+S (B) a 


2.97 


free 


32.41 


4.82 






34.5 


0.92 


11.55 


-20.91 


3.9 


4.4 


0.83 


1.12 


S+S (D) a 


1.90 


free 


2.38 


0.35 






29.4 


0.88 


13.31 


-19.16 










NGC 4694 S 


3.53 


free 


32.92 


2.52 






-34.9 0.46 


11.23 


-19.78 


2.6 


3.6 




1.20 


[69.8,187.5] S+E(B) 


3.23 


free 


25.84 


1.98 






-34.3 0.42 


11.49 


-19.52 


2.8 


4.3 


0.82 


1.06 


S+E (D) 


1.00 


fix 


51.23 


3.93 


30.53 


2.34 


-62.5 0.79 


13.12 


-17.89 










DeV 


4.00 


fix 


37.93 


2.91 






-34.9 0.46 


11.17 


-19.84 


1.7 


1.6 




1.29 


S+S (B) 


3.25 


free 


27.12 


2.08 






-34.5 


0.43 


11.42 


-19.59 


3.1 


4.4 


0.87 


1.05 


S+S (D) 


0.61 


free 


52.81 


4.05 






-82.6 0.79 


13.51 


-17.50 










NGC 5018 S 


5.39 


free 


23.83 


4.41 






-83.4 


0.71 


9.99 


-22.96 


2.6 


2.6 




1.26 


[19.6,247.5] S+E(B) 


5.29 


free 


19.08 


3.53 






-84.4 


0.67 


10.22 


-22.73 


4.3 


4.1 


0.87 


1.08 


S+E (D) 


1.00 


fix 


30.73 


5.69 


18.31 


3.39 


-30.8 


0.90 


12.26 


-20.68 










DeV 


4.00 


fix 


16.93 


3.13 






-84.1 


0.72 


10.12 


-22.83 


8.9 


8.2 




2.20 


S+S (B) 


5.13 


free 


11.23 


2.08 






-86.6 


0.64 


10.82 


-22.13 


4.7 


4.5 


0.51 


1.06 


S+S (D) 


2.45 


free 


31.55 


5.83 






-75.8 


0.78 


10.86 


-22.09 










NGC 5576 S 


6.11 


free 


31.07 


3.05 






87.7 


0.70 


10.10 


-21.45 


4.9 


4.9 




1.48 


[11.8,247.5] S+E(B) 


3.93 


free 


11.82 


1.16 






89.2 


0.69 


10.68 


-20.88 


4.7 


5.1 


0.65 


1.05 


S+E (D) 


1.00 


fix 


78.84 


7.74 


46.98 


4.61 


79.9 


0.72 


11.33 


-20.22 










DeV 


4.00 


fix 


17.44 


1.71 






87.9 


0.72 


10.33 


-21.23 


17.1 


17.0 




3.61 


S+S (B) 


3.45 


free 


8.89 


0.87 






89.6 


0.69 


10.90 


-20.66 


4.4 


4.7 


0.52 


1.04 


S+S (D) 


1.48 


free 


71.10 


6.98 






82.2 


0.72 


10.98 


-20.57 











NOTE. — Col. (1): Object name followed by inner and outer boundaries in arcsecs within square brackets, used to estimate the luminosity of 
the tidal features. Col. (2): The models used to fit galaxies with GALFIT. S represents the Sersic, E represents the exponential and DeV is the de 
Vaucouleurs profile. B and D in parentheses denote the bulge and disk component, respectively. Col. (3): Sersic index. Col. (4): Type of the fit while 
running GALFIT for each component. The Exponential and de Vaucouleurs profile are fixed to have n=l and 4, respectively and denoted as "fix". 
Col. (5): Effective radius measured along the semi major axis in arcsecs for bulge and also for disk. Note that for exponential disks, the scale length 
can be obtained by dividing this effective radius by 1.678. Col. (6): Effective radius measured along the semi major axis in kpc, estimated with the 
heliocentric redshift listed in NED. Col. (7): Disk scale length in arcsecs for disk component. Col. (8): Disk scale length in kpc for disk component. 
Col. (9): Position Angle, east of north. (North = 0, East = 90) Col. (10): Minor axis to major axis ratio. Col. (1 1): Apparent AB magnitude at 3.6 
/^m. Col. (12): Absolute AB magnitude at 3.6 [im. Col. (13): The ratio of the residual to the galaxy model in the circular region between inner and 
outer boundary. Note that in the case of Sersic with exponential fit, we estimate the ratio with the total model of bulge and disk, but we put the value 
just in the first row in the S+E section. Col. (14): The ratio of the residual to the galaxy model in an elliptical boundary using PA, b/a, inner and outer 
boundary. Col. (15): Bulge to Total ratio. Col. (16): The reduced \ 2 , which is a indicator of goodness of fit. 

a From the Sersic+Sersic, two-component fit, we took the larger Sersic index component as a bulge component. But in this case (NGC 2634), the 
effective radius of the larger Sersic index component (bulge) is larger than that of smaller Sersic index component (disk), hence this model fit may 
not be realistic. 



effective radius and the mean surface brightness at that radius 
for ETGs showing tidal features in Figure 9. To make compar- 
isons with other studies, we also plot KRs from the literature, 
which were obtained from the same 3.6 /im-band. Each panel 
shows the KR obtained with a specific model fit. We show the 
KR derived from the n=4 profile fit in the uppermost panel to- 
gether with the results of ETGs from Jun & Im (2008) in grey 
triangles and Falcon-Barroso et al. (201 1) in red squares. Note 
that structural parameters (r e ff and < /i e ff >) of these two stud- 
ies were obtained with one-dimensional R l l A profile fittings, 
while our data were obtained with two-dimensional model 



fits. ETGs presented in Falcon-Barroso et al. (201 1) can be di- 
vided into two groups that have Sersic index n=4 and n=free. 
We only plot ETGs that have n=4 in this panel. In the second 
panel of Figure 9, we plot the result from a Sersic (n=free) 
fit. To compare with normal elliptical galaxies, we draw a red 
dotted line adopted from Figure 8 of Fisher & Drory (2010), 
which represents elliptical galaxies of the Virgo cluster from 
(Kormendy et al. 2009). Since these galaxies were originally 
studied in the optical band, Fisher & Drory (2010) shifted the 
surface brightness measurement to the L'-band. We also plot 
ETGs from Falcon-Barroso et al. (201 1) that were fit with Ser- 
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FIG. 8. — Comparison of structural parameters among different choices of models. The numbers next to each point are the NGC numbers of the sample galaxies 
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correspondence. Note that the Sersic+Sersic model fit of NGC 1222 and NGC 2634 may not be realistic, as denoted in Table 4. 



sic index n=free in orange squares. In the third panel, we show 
the result of a Sersic (for the bulge) plus exponential fit (for 
the disk). In the Sersic+exponential fit, the values of r et£ and 
/x e ff are those of the bulge component. We also plot the clas- 
sical bulge (in orange triangle) and pseudo bulge (in green 
triangle) of Fisher & Drory (2010). Their structural parame- 
ters were determined by fitting the surface brightness profiles 
with a one-dimensional Sersic function plus an exponential 
outer disk. The last panel shows the result of the fit with a 
Sersic plus Sersic function. The classical and pseudo bulges 
are the same as above, and they were plotted just for com- 
parison. However, because a different fitting model is used, 
we can not apply direct comparison in the lower most panel. 
NGC 4694 is a galaxy with a nuclear HII region (Veron-Cetty 
& Veron 2006) and the velocity dispersion of this galaxy is 61 
km s" 1 . Therefore it falls below the main KR of ETGs, simi- 
lar to dwarf galaxies showing nuclear star formation. Except 
for NGC 4694, ETGs with tidal features do not stand out in 
the KR from the non-interacting normal ETGs of Jun & Im 
(2008) and Falcon-Barroso et al. (201 1). 

We present the FR which is an empirical scaling relation 
between the effective radius, the mean surface brightness, and 
the central velocity dispersion of tidally disturbed ETGs in 
Figure 10, where the axes are projected in the direction of the 
smallest scatter at 3.6 \im as provided by Jun & Im (2008). 
As in Figure 9, we show FPs derived with four different sur- 
face brightness profile models. In this figure, we compare 
the FPs of our sample galaxies only with those of Jun & Im 
(2008). The triangles of Jun & Im (2008) represent the re- 
sults obtained with the R X I A profile, hence only the uppermost 
panel shows a direct comparison, but we plot them together 
in the other panels for comparison. Velocity dispersions were 
adopted from Hyperleda. As discussed in the previous sec- 
tion, the Sersic+Sersic model fit of NGC 1222 and NGC 2634 
may not be ideal because the effective radii of the bulge com- 
ponent are larger than those of the disk component. So, we 
draw grey arrows in the lowermost panel of Figure 9 and 10 
representing an exchange of the values of r e ff and /i e ff of the 
bulge and disk components. We find that tidally disturbed 
ETGs lie on the same plane as do normal ETGs (Jun & Im 
2008) and Virgo ellipticals (Kormendy et al. 2009). Because 
galaxies exhibiting tidal features do not occupy a distinct re- 



gion in the KR or FR this implies that the structural properties 
of the spheroidal components of tidally disturbed galaxies are 
not significantly different from those of the other ellipticals or 
the spheroids of galaxies with no signs of tidal interaction. 

Our result is consistent with the findings of numerical sim- 
ulations of mergers that show that the FP is preserved dur- 
ing gas-poor, dissipationless major merging events (Capelato 
et al. 1995; Gonzalez-Garcia & van Albada 2003; Nipoti 
et al. 2003; Boylan-Kolchin et al. 2005; Boylan-Kolchin et al. 
2006). Rothberg & Joseph (2006) also find observationally 
that most of the single-nuclei merger remnants, including 
shell galaxies, appear to lie on the TT-band FP or within the 
scatter of the ellipticals taken from Pahre (1999). 

On the other hand, simulations show that projections of 
the FP, such as the KR and the Faber-Jackson relation (Faber 
& Jackson 1976, FJR), are not preserved via dry mergings 
(Nipoti et al. 2003; Boylan-Kolchin et al. 2006). Boylan- 
Kolchin et al. (2006) showed that the location of merger rem- 
nants in the FJR and KR is sensitive to the amount of the 
energy transferred from the bulge to the dark matter haloes 
during the merger. They claimed that mergers on orbit with 
large angular momentum would experience more dynamical 
friction. This results in galaxies with larger velocity disper- 
sion (smaller f3 in /? e ff °c Lcr^) and smaller r e ff (smaller a in 
R e ff oc M"). Therefore, such galaxies form a less steep KR 
and FJR. While the end product of mergers on orbits with less 
angular momentum will form steeper KR and FJR, i.e., have 
higher a and (5. But in Figure 5 (KR) of Rothberg & Joseph 
(2006), shell galaxies and other merger remnants, except for 
luminous and ultra-luminous infrared galaxy mergers, also lie 
within the scatter (Pahre 1999), consistent with our result in 
the KR. 

Our results do not mean that the properties of the galax- 
ies showing tidal debris have not changed as a result of the 
merging events. It is more likely that tidal features induced 
in minor mergers or interactions, do not affect bulge proper- 
ties. Alternatively if those tidal features were caused by major 
mergers, then enough time has passed for their bulges to relax 
and become similar to those of ETGs that do not show tidal 
disturbances. By checking H/?o line indices, Falcon-Barroso 
et al. (201 1) show that outliers from the FP tend to have young 
stellar populations in common. Having young stellar popula- 
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tions might be a more efficient way to push galaxies off the 
FP than dry merging as also suggested by the E+A findings 
(Yang et al. 2008). We infer that the dynamical scattering of 
the FP due to a merging event last only for a relatively short 
time scale. This is also in agreement with the finding that 
ETGs show a tight relation in the FP. Thus through dry merg- 
ing events between galaxies which do not have enough gas to 
trigger star formation, galaxies appear to slide along the FP in 
the direction of larger effective radii, rather than lie off the FP. 

6. SUMMARY AND CONCLUSIONS 



In a sample of ETGs drawn from S 4 G that have T < 0, 
we perform unsharp masking, two-dimensional fitting, and 
derive structure maps to search for tidally disturbed ETGs. 
We identify shells, ripples, tidal tails or broad fans of stars 
in 17 ± 3% (1 1/65) of ETGs with the imaging data that have 
average depth of < fj,%64B > (lc) ~ 27 mag arcsec -2 at 3.6 
fjm. 

We estimate that the tidal disturbances contain 3 - 10% 
of the total 3.6 pm galaxy luminosity, which can be directly 
translated to the stellar mass contained in the debris. 

To investigate the structural properties of tidally disturbed 
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galaxies, we modeled the galaxies with four different func- 
tions, de Vaucouleurs, Sersic, Sersic+exponential and Ser- 
sic+Sersic profile. We find that the precise measurement of 
structural parameters depends on the adopted decomposition 
model. We also find that the total magnitudes of galaxies vary 
among different model fits. 

We explored the position of these galaxies in the Kormendy 
relation and Fundamental Plane and find that ETGs with tidal 
debris occupy the same area as those galaxies without tidal 
debris or bulges of early type systems. This implies that the 
bulge properties of tidally disturbed galaxies are not signifi- 



cantly different from those of normal ETGs. We infer that dry 
merging cause galaxies to slide along the FP and minor merg- 
ing or weak interactions that cause tidal disturbances in our 
sample either affect only outer part of the galaxy, or enough 
time has elapsed that the bulges of the tidally disturbed galax- 
ies are well relaxed. Thus, having young stellar populations 
might be a more efficient way to make galaxies deviate from 
the FP than going through dry merging events. 

The S 4 G will survey over 2,300 nearby galaxies by the end 
of 2012, and this study will be extended to include a complete 
set of 180 ETGs. 
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